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We study the coupled dynamics of spin and charge currents in a two-dimensional electron gas in 
the transport diffusive regime. For systems with inversion symmetry there are established relations 
between the spin Hall effect, the anomalous Ifall effect and the inverse spin ffall effect. However, in 
two-dimensional electron gases of semiconductors like GaAs, inversion symmetry is broken so that 
the standard arguments do not apply. We demonstrate that in the presence of a Rashba type of 
spin-orbit coupling (broken structural inversion symmetry) the anomalous Hall effect, the spin Hall 
and inverse spin Hall effect are substantially different effects. Furthermore we discuss the inverse 
spin Hall effect for a two-dimensional electron gas with Rashba and Dresselhaus spin-orbit coupling; 
our results agree with a recent experiment. 



Despite the anomalous and spin Hall effect being 
closely related, their histories are rather different. The 
anomalous Hall effect was experimentally discovered [T] 
almost at the same time as the ordinary Hall effect, while 
the spin Hall effect, first predicted in 1971 [21 13] and sev- 
eral times recently [1H7], has been experimentally seen 
only in the last few years f51111|. This is not surprising 
as the anomalous Hall effect entails the measurement of 
currents and voltages which is well established experi- 
mentally, whereas the spin Hall effect requires the detec- 
tion of a spin current, which has to be done in an indirect 
way; for a review see for example [TH [T3] and [HI [T5] . 

Since the anomalous and spin Hall effect have the same 
physical origin, namely the spin-orbit interaction which 
couples charge and spin degrees of freedom, their de- 
pendence on various physical parameters is expected to 
share similar trends. Depending on whether the spin- 
orbit coupling is intrinsic in the band structure or ap- 
pears due to coupling to impurities one speaks about 
intrinsic or extrinsic mechanisms. The interplay of in- 
trinsic and extrinsic mechanisms is nontrivial. For in- 
stance the intrinsic Rashba type of spin-orbit coupling in 
a two-dimensional electron gas suppresses drastically the 
extrinsic (skew-scattering) contribution to the spin Hall 
conductivity [TMI^ . 

It is the purpose of this paper to develop a similar 
analysis for the anomalous Hall effect and for the inverse 
spin Hall effect. We will start with a phenomenological 
discussion of charge and spin currents in a metal or semi- 
conductor with diffusive charge carrier dynamics. We will 
study in detail the two-dimensional Rashba model includ- 
ing extrinsic skew-scattering. In the anomalous Hall con- 
ductivity we find an unexpected anomaly in the magnetic 
field dependence. Our analysis of the inverse spin Hall 
effect for a system with Rashba and Dresselhaus spin- 
orbit coupling is consistent with the experimental results 

of [20]. 

As a starting point we consider a system with spin- 



orbit coupling, where a spin polarized current in the x- 
direction generates a small current 6jy into the transverse 
direction with 

Clearly, from these equations we can conclude that: a) a 
spin polarized current generates a transverse charge cur- 
rent (anomalous Hall effect), b) a pure charge current 
{jxi = jxi) generates a transverse spin current (spin Hall 
effect), and c) a pure spin current [j^i; = —jxi) gener- 
ates a transverse charge current (inverse spin Hall effect). 
The three effects mentioned arc thus very closely related, 
and the magnitude of all of them is determined by the 
dimensionless parameter 7. 

Often, however, Eq. ([!]) is not sufficient for the theoret- 
ical description and, in the following, we will achieve the 
necessary generalization of the equations. Let us write 
the current in the x-direction as 

jx-\ = o-^Ex^, = (TiEx\^, (2) 

where a^^i and Ex^^i are the spin-dependent conduc- 
tivity and electric field in the x-direction, respectively. 
In order to allow later arbitrary directions of the spin- 
polarization, we find it convenient to introduce here the 
charge and spin components for the field and the current, 
Ex^,i ^ Ex± \Exz and j^^^^ = \jx ± jxz- Equation ^ 
can now be rewritten as 

jx = O-Ex + (TOzExz, (3) 
jxz = '^(^Exz+ (TzqEx, (4) 

where a — cr^ + = /ip is the Drude conductivity, p 
and /i being the charge density and the mobility, re- 
spectively. The conductivity oqz = IJ'Sz, with Sz the 
spin density, mixes spin and charge currents and ap- 
pears due to the fact that electrons carry both degrees 
of freedom; the Onsager relations require the symmetry 
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o'oz(sz) = —crzo{—Sz)- Notice that the charge and spin 
currents (as weU as charge and spin density) as defined 
here have equal units. The transverse currents are given 
by 



+ lOzJx, 
lOzjxz, 



(5) 
(6) 



with 7 = |(7t + 7^) E^nd = (7-t- ~ 7^) when we aUow 
different 7s for spin up an down. 

In the next step diffusive currents are considered too. 
This is achieved be replacing the electric fields by 



tCtEx 
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<j£xz = -DdxSz 



:<7Exz 



(7) 

(8) 



The diffusion coefficient, D, is related to the conductiv- 
ity via the relation a = 2e^DNQ, where Nq is the single- 
particle density of states at the Fermi energy. Allowing 
now an arbitrary direction of the fields and the spin po- 
larization we obtain the set of equations 



Ji 

Jla 
Sjla 



(9) 

^lablabia, (10) 



1 



= IQab [f^h + croc^bc] + ^lablOcjbc 



(11) 

(12) 



The structure of Eqs. ([9]) - ( 12 1 is similar to the equations 
given in Ref. [H], the difference being the terms a^a^ia 
and 7o6 ja in the charge current and jocjbc in the spin one, 
which do not appear in [5T]. The last term in Eq. (12) 



will however be of no importance in the present article 
and, as such, will be ignored in the following. 

We proceed by calculating the parameters entering 
Eqs. ([9]) - (fnl) from a microscopic model. We consider 



a disordered two-dimensional electron gas (2DEG) with 
Hamiltonian 



H 



P 

2m 



IK 
h 4 



1. 



pcr X axV(x) p -h (T. 

(13) 

In this Hamiltonian we have both intrinsic and extrinsic 
spin-orbit coupling. The intrinsic spin-orbit interaction 
modifies the band-structure and enters in the form of a 
spin-dependent vector potential [T^[2M25] . which for the 
Rashba model is given by 



ma 

A — —— (T X Bz 

h 



(14) 



with the only non-zero components A^y = —Ayx = 
2ma/h. V(x) is the scalar potential due to the scattering 
from impurities and gives rise to the extrinsic spin-orbit 
coupling with strength characterized by the length Aq. 
Both spin-orbit couplings are assumed to be weak, i.e. 



A • pi?/m ^ and XqPf ^ ^- The Zeeman field b may 
be due to an external magnetic field or may arise due to 
the exchange field of a ferromagnet. In the following, for 
the sake of simplicity, we take units such that h — \. 

For our microscopic model the density of states is 
A^o = m/2-K and the diffusion constant is D ~ ^VpT, with 
T the elastic scattering time. The latter is determined 
from the disorder potential and in the Born approxima- 
tion, assuming (y(x)V(x')) = (5(x - x')/(27rAror). The 
parameter 7 has, in principle, contributions from the 
skew-scattering, side-jump, and the intrinsic mechanism. 
In this paper, motivated by the fact that in 2DEGs the 
skew-scattering is typically considerably stronger than 
the side-jump, we limit our discussion to the interplay 
of skew-scattering and the intrinsic mechanism. We then 
write the parameter 7 as 



7 = 7skow + 7int , 



with 



7skcw 
7int 



16 
-ma^T. 



(15) 

(16) 
(17) 



For an explicit derivation, one may see [ini [2S]- In 
Eq. ( 16 ) vq is the scattering amplitude from the im- 
purity potential. For the parameter 7oz we find in 
our model only a skew-scattering contribution, explic- 
itly jqz — 47si5cwfoz/f ■ It is stressed in Ref. [21^ that 
Eqs. ^ - ( 12 1 are only valid in systems with inversion 
center. In the absence of the inversion symmetry - which 
is the case in the situation we consider here - extra terms 
appear. However for our model Hamiltonian ( 13 1 these 



extra terms are conveniently taken into account by a re- 
definition of the field £ia, which is now given by 
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aSia = -Ddiisa - si'') - DeabcAlbis, - s^«), (18) 



where s*^' = (— e)^Aob. For example the fields £yz and 
£xz are given by 



:<y£x 



-Ddy{sz 
-Ddx{sz 



s:^)-fi?(2ma)(5,-4«),(19) 
s7) + D{2ma)(sx-s7).{2Q) 



Again we refer to the literature for microscopic deriva- 
tions. For example in |27j the expressions for the spin 
and charge currents in the case Aq = were obtained by 
exploiting an SU(2) symmetry of the Rashba model. 

How does our approach compare with other studies 
of the diffusive dynamics of spin and charge? Combin- 
ing the current density ^ and (111 with the continuity 
equations for spin and charge one finds coupled diffusion 
equations. Such diffusion equations have been derived 
for the Rashba model e.g. in [551 US] and for the system 
with both a Rashba and a linear Dresselhaus term in [3D] . 
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Our analysis extends these works: whereas the cited pa- 
pers concentrate on the intrinsic spin-orbit couphng we 
include also the experimentally relevant skew-scattering. 
Furthermore the spin-charge coupling conductivities croa 
and (Tao are neglected in pSHSO] . In the following we will 
apply the formalism to the various Hall effects. 

Anomalous Hall effect and spin Hall effect: The 
anomalous Hall effect describes a contribution to the 
Hall conductivity due to the spontaneous magnetization 
in a ferromagnet, the Hall current being perpendicular 
to both the magnetization and the electric field. The 
spin Hall effect consists instead in the appearance of a 
spin current orthogonal to an applied electric field in 
a non-magnetic material. Let us assume homogeneous 
conditions, take the electric field along the x-axis and 
the magnetization along the z-axis, and write down the 
charge and spin currents along the y-axis. To linear order 
in the electric field Eqs. ^ and ( |TT| ) become 

jy = (70z£yz + + -I- 7skew)o-zO-E2: (21) 

jyz = '^<j£yz+lCFEj:+"iao:,£xz- (22) 

Let us at first examine the anomalous Hall current, 
Eq. (21 1, in the pure Rashba model (Aq = 0). It is known 



that in the presence of spin-orbit coupling an electric field 
induces a spin polarization [31 . In the Rashba model this 
lies in-planc, and for our geometry along y [32 , which 
implies that E^z — 0, but £yz ^ . From Eqs. (17) and 
( 19 1 we get 



4D(2r 



-Sy — Ama tE^ 



(23) 



For the spin polarization one has 



Sy = e^NoarE^, 



and thus the anomalous Hall effect in the pure Rashba 
model vanishes in agreement with explicit diagrammatic 
calculations. Notice that in the diagrammatic calcula- 
tions a finite anomalous Hall effect is found from a skew- 
scattering-like contribution which appears in a higher or- 
der in the magnetic field and in the presence of magnetic 
impurities |33H36j . Such a contribution we do not con- 
sider here. 

The disappearance of the anomalous Hall effect is re- 
lated to the vanishing of the spin Hall effect in the pure 
Rashba model. To see this, let us consider the spin Hall 
current (22 1. Since = 0; we find by comparing Eqs. 



(211 and (22) the relation 



Jjf — * ^ Jyz, 



(25) 



so that a vanishing spin Hall current implies a vanishing 
charge Hall current. 

This relation is no longer true in the presence of 
both intrinsic {a ^ 0) and extrinsic spin-orbit coupling 



(Ao 0). In this case, the combination of an out-of-plane 
magnetic field or exchange field together with an in-plane 
electric field (in the x-direction) generates a component 
of the spin-polarization in the x-direction so that the field 
Sxz no longer vanishes. 

To calculate the spin polarization we borrow from |19j 
the equations 



f(s 



Gq J 



boff X s + Sb. 



(26) 



Here f is the spin relaxation matrix which in the case of 
pure Dyakonov-Perel spin relaxation reads 

— diag(l,l,2), l/TDP^D{2ma)\ (27) 
tdp 

The spins relax towards the equilibrium density Sg^ = 
— e^A'pb and precess in an effective magnetic field. 



Jeff 



b + 2maiiez x E, 



(28) 



whereas is an electric field dependent source term 



= 2TOacr(7skow + 7int)ez x E. 



(29) 



Solving these equations in the static limit and ignoring 
possible nonlinearities in the electric field, the spin po- 
larization is determined as 



&zT^ 



1 + {hzTDpf 



(30) 



e NqutE^ - , , ^7skew-e;a; (31) 



1 + {hzTDP? 



(32) 



Knowing the spin polarization we can now calculate the 



(24) Hall and spin Hall current, using Eqs. (21) and (22). In 



the weak magnetic field limit [b^Top ^ 1) the result for 
the Hall current is 

jy= ( 1 + i^^iS^ ) %-1s]icv,(yOzEx, (33) 
V 2 7iiit / 

which means that a weak Rashba term (7skew ^ 7int) 
may considerably enhance the anomalous Hall effect. 
The spin Hall current, however, vanishes 157] 



Jyz 







(34) 



in the presence of the Rashba coupling. The term 



7skGw/7int on the right hand side of Eq. (33) appears to 
be singular when the Rashba coupling goes to zero. This 



is because we have assumed in Eq. (27 1 that the Elliott 



Yafet spin relaxation rate, I/ts, can be neglected when 
compared to the Dyakonov-Perel one, i.e., Ts 3> Top. 



The same holds in Eq. (34). Clearly this assumption 



is no longer valid when a goes to zero as it has been 
discussed in Ref. [12] in connection with the spin Hall ef- 
fect. In order to include the Elliott- Yafet relaxation due 
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FIG. 1: The Hall angle as a function of the magnetic field. 
Since we concentrate on the anomalous Hall effect the con- 
tribution due to the Lorentz force is ignored. We estimate 
the Hall angle using parameters valid for GaAs with a car- 
rier density of 10^^/cm'^. In the absence of Rashba spin-orbit 
coupling the Hall angle is determined from skew-scattering, 
with 7skew ~ 2.7 X 10~^ assuming positively charged impu- 
rities {NqVo < 0). We obtained the full line assuming a mo- 
bility of /i = 10*cm^/Vs and a spin-orbit coupling constant 
a = 10~^'^eVm. The long-dashed and dashed lines correspond 



to /I = 10 cm /Vs, a 
a = 10""eVm. 



10-"eVm and ^ = IQ-^cmVYs, 



to the extrinsic spin-orbit interaction, we have to modify 



Eq. ( 27 ) in the foUovifing way [HI 
^ 1 



1 



diag(l,l,2) + -diag(l,l,0). (35) 
tdp Ts 



As a consequence, in Eqs. ( 30]|3T ) we must operate the 
replacement r^p — )■ r^p -I- r^T^^ which guarantees the 
correct a limit. 

In the strong magnetic field limit {hzTop ^ 1), and 
again assuming ^ tdp, the anomalous Hall current is 
given by 



jy — SjskcwCTQzEx 



(36) 



which is identical to the result in the absence of Rashba 
spin-orbit coupling. 

The Hall angle, jy/jx, as function of the magnetic field 
is shown in Fig. [T] for different values of the mobility 
and the Rashba term. In the absence of Rashba spin- 
orbit coupling Eq. ( [36| ) implies that jy/jx = ^jskcwb/sp, 
i.e. the Hall angle as a function of the magnetic field is a 
structureless line. The Rashba term causes an anomaly 
in weak magnetic fields. The width of this anomaly is 
set by the Dyakonov-Perel relaxation rate, and therefore 
depends strongly on the value of the Rashba coupling but 
also on the mobility. 

Inverse spin Hall effect. In the inverse spin Hall effect 
(ISHE) spin-polarized carriers are injected into a non- 
magnetic material. In Ref. [10] . for instance, the injec- 
tion of the spin current was achieved via a ferromagnet 
contacting the spin-orbit active material, while in Ref. 
[20j spins were injected by applying an optical technique. 



The spin-current generates a transverse charge current 
which in the end is detected via a standard voltage mea- 
surement. We consider the situation where a spin current 
is injected in the x-direction and generates a charge cur- 
rent in the y-direction. We analyze the ISHE via Eq. ^ 
assuming jy is linear in the driving force. Then the ex- 
pression for the current considerably simplifies and reads 



jy = '^ijx 



(37) 



From this relation one can directly read off a Hall an- 
gle. If a fully polarized current is injected into the system, 
then jxz = \ jx at the injection point and the Hall angle 
is given by 



"if ~ jy/jx = 27- 



(38) 



For a — 0, i.e. without the intrinsic spin- 
orbit coupling, the spin density decays exponentially 
with the distance from the injection point, Sz(x) — 
Sz(0) exp(— x/ig), the spin current is proportional to the 
derivative of the spin density, jxz = ^DdxSz, and there- 
fore the Hall angle drops exponentially too. In [501 where 
both intrinsic and extrinsic spin-orbit coupling is present, 
a spin profile of the type Sz{x) w Sz{0) cos{Qx) is ex- 
pected, where the constant Q depends on the strength of 
the intrinsic spin-orbit coupling. Surprisingly the mea- 
sured Hall data is consistent with the assumption of 
jy being proportional to the spin density instead of its 
derivative. In the following we analyze the experiment in 
more detail in order to understand this point. The ex- 
periment was designed such that the linear Rashba and 
Dresselhaus spin-orbit coupling in the 2DEG are of sim- 
ilar size. For simplicity we assume that both terms are 
equal, so that the Hamiltonian is given by (a — — /3) 



H : 



2m 



+ a{py -Px){(^x + (^y)- 



(39) 



It is useful to formulate the theory in a rotated frame with 
unit vectors e-|_ = (ej;+ey)/%/2 and e_ ~ (ej.— ej,)/%/2 so 
that the spin-orbit coupling reads a{py —Px)i<^x '^'^y) — 
— 2Q:p_cr_|_. Solving the spin diffusion equation with the 
boundary condition s(0) — Sz{0)ez we find a spin-spiral 
of the form 




.(O)e' 





— sin(4mQ:a;- 
cos(4mQ;x_ 



(40) 



the persistent spin helix of [30l [38]. Here the spin- 
relaxation length Lg was introduced by hand but can 
be justified microscopically by any spin-relaxation mech- 
anism like e.g. the Elliott- Yafet mechanism. In the latter 
case, one finds explicitly Ls = \/2Dts if 3> t/^p. In 
[20] the 2DEG channel is patterned along the [110] di- 
rection, the direction of the spin-helix propagation. The 
Hall current is then proportional to the spin current flow- 
ing in the a;_-direction. After modifying Eq. ( 19 1 to 
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include both a Rashba and a Dresselhaus term (with 
a = — /3) we find 

jy = 47j^_z =47[-£»a^_s^ + £)(4TOa)s_] (41) 

= 47-^5^(0) cos(4TOaa;_)exp(-a;_/Ls). (42) 

The Hall current indeed follows the spin polarization in 
the 2:-direction with periodic changes of the sign with in- 
creasing distance from the spin-injection point, in agree- 
ment with the experimental finding. Also the absolute 
value of the Hall angle, which is of the order of some 
10~^ is consistent with realistic estimates of the param- 
eters. 

Summary: We presented equations describing the cou- 
pled dynamics of spin and charge currents in a two- 
dimensional electron gas. Unlike in inversion symmetric 
systems, where the spin Hall effect, the anomalous Hall 
effect and the inverse spin Hall effect are essentially the 
same thing no general relation between the three effects 
can be given. For example in the pure Rashba model we 
find a vanishing spin Hall and anomalous Hall effect, but 
a finite inverse spin Hall effect. We analyzed the inverse 
spin Hall effect for a system where Rashba and Dressel- 
haus spin-orbit coupling have equal strength; our results 
compare well with a recent experiment. 

This work was supported by the Deutsche Forschungs- 
gemeinschaft through SPP 1285 and partially supported 
by EU through PITN-GA-2009-234970. 
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